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Kazuto Noda0 Robert Peters, and Norio Kawakami 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 



(N 

o 

(N 



oo 



CZ2 

wo; 

•4-* 

o3 



c 

o 
o 



> 

oo 
in 

in 
o 

(N 



X 



Thomas Pruschke 

Institut fur Theoretische Physik, Universitat Gottingen, Gottingen D-37077, Germany 

We investigate many-body effects on a mixture of interacting bosons and fermions loaded in 
an optical lattice using a generalized dynamical mean field theory combined with the numerical 
renormalization group. We show that strong correlation effects emerge in the presence of bosonic 
superfluidity, leading to a renormalized peak structure near the Fermi level in the density of states 
for fermions. Remarkably, this kind of strong renormalization appears not only in the metallic phase 
but also in the insulating phases of fermions such as in the empty/filled band limit. A systematic 
analysis of the relation between the quasiparticle weight and the strength of superfluidity reveals 
that the renormalization effect is indeed caused by the boson degrees of freedom. It is found that 
such renormalization is also relevant to a supersolid phase consisting of a density wave ordering of 
fermions accompanied by bosonic superfluidity. This sheds light on the origin of the peak structure 
in the supersolid phase. 

PACS numbers: 



I. INTRODUCTION 

Strongly interacting atoms in optical lattices have at- 
tracted much interest [![. In these systems, one can tunc 
the interaction strength and the lattice structure by con- 
trolling a magnetic field and the intensity of lasers. Due 
to such high controllability, the cold atom systems can be 
regarded as efficient simulators of quantum many-body 
physics. Several intriguing phenomena related to strong 
correlations were already observed experimentally, for ex- 
ample, a superfluid-Mott transition in bosonic systems 
, a metal- Mott transition in fermionic systems [H, ||j , 
etc. 

Experimental research has already been extended to 
the topics which are not easy to investigate in conven- 
tional condensed matter physics. A typical example is a 
Bose-Fermi (BF) mixture realized in a harmonic trap 0- 
[Toj . and also in optical lattices [Tll - [l7j . Rapid progress 
in these experiments has stimulated theoretical research 
on related topics [l8l - [27| ; for example, the possibility of 
fascinating ground states such as a supersolid state (co- 
existence of bosonic superfluidity and densit y-w ave or- 
dering of fermions) has been proposed 18 L 22^2414261. In 
this context, Titvinidze et al. in Refs. [24, l25l| pointed 
out that the density of states for fermions has an anoma- 
lous peak near the Fermi surface in the supersolid phase. 
It should be important and interesting to further clarify 
the many-body aspects of this structure, which naturally 
motivates us to provide a detailed analysis of it. 

On the other hand, dynamical properties of the BF 
mixture systems have not been studied systematically. 
Recent rapid advances in probing dynamical properties 
of fermions via the rf spectroscopy make it possible to 
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observe the single-particle excitation spectrum in a BCS- 
BEC crossover region [28 30], the quasiparticle weight in 
a polaronic system (an extremely imbalance system) [3ll . 
|32j etc. These experimental developments would be also 
applied for the BF mixture systems in the near future, 
encouraging theoretical studies on dynamical properties 
of the BF mixture systems. 



Motivated by these research activities, we investigate 
many-body effects on the BF Hubbard Hamiltonian with 
particular emphasis on its dynamical properties. We re- 
veal unique features inherent in the BF mixture systems; 
the renormalization effect caused by the boson degrees of 
freedom gives rise to a peak structure near the Fermi level 
for the density of states both for metallic and insulating 
phases. A similar peak structure due to many-body ef- 
fects appears even in the supersolid phase, which we will 
discuss in the following sections. Our systematic study 
in this paper clearly explains that the origin of the peak 
structure is many-body effects induced by the interplay 
of the boson and fcrmion degrees of freedom. 

This paper is organized as follows. In Sec. HH we 
introduce the BF Hubbard model, and briefly explain 
the framework of a generalized version of dynamical 
mean field theory (DMFT), which extends the original 
fermionic DMFT to treat additional boson degrees of 
freedom. We make use of the numerical renormalization 
group (NRG) as an impurity solver of DMFT. In Sec. [TTI1 
we reveal the renormalization effects in the presence of 
bosonic superfluidity In the last part of this section, we 
shortly discuss how the many-body effects occur in the 
impurity model corresponding to the effective impurity 
model in the DMFT calculation. In Sec lIVI we discuss 
how the renormalization effects appear in a supersolid 
state. In Sec. [Vj we summarize our results. 
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II. MODEL AND METHOD 

We consider a BF mixed system, which can be de- 
scribed by the following BF Hubbard Hamiltonian, 



H = H b + H f +H bf , 



(1) 



(i,3) 



f n{, 



(2) 
(3) 

(4) 



with n b = b\bi (n( = /,•/,), where &i(/,) annihilates a bo- 
son (fermion) at site i. Here, t b (t^) is the boson (fcrmion) 
transfer integral, /i b (^) the chemical potential for bosons 
(fermions) and U b (U b *) the on-site boson-boson (boson- 
fermion) interaction. Note that denotes the summa- 
tion over the neighboring lattice sites. Note that we treat 
spinlcss (one-species) fermions in this article. The cor- 
responding circumstance, where one species of fermions 
are mixed with bosons, has already been realized experi- 
mentally by properly selecting one of the hypcrfinc states 
in fermionic atoms (see @). 

To investigate the ground-state as well as dynamical 
properties of the system, we employ the DMFT. In or- 
der to treat the boson degrees of freedom, we use a gen- 
eralized version of DMFT which is introduced in Ref. 
[HHl]. In the generalized DMFT, the lattice model is 
mapped onto an effective impurity model embedded in 
an effective medium, as usually done in DMFT [HI, HBJ . 
The Green's function is obtained via the self-consistent 
solution of this impurity model. This is why the DMFT 
exactly includes local quantum fluctuations, which can- 
not be taken into account by conventional mean-field ap- 
proaches. The different point from the fermionic DMFT 
is that the bosonic superfiuid order parameter, ip — (b), 
should be obtained in self-consistency steps. We perform 
the calculation using a semielliptic local density of states 
with a bandwidth W = 4i* for the noninteracting sys- 
tem. In the following, the half bandwidth D = 2t' is 
used as a unit of energy. 

In order to discuss many-body effects in the phases 
with and without a density-wave order, we introduce the 
corresponding effective impurity models. When we ana- 
lyze the phases without a density wave order, we use the 
following generalized single impurity Anderson model as 
an effective impurity model, 



6 b 

n 



ffcsiAM = -zt\ytf + p*b) + ^-n b (n b - 1) - A t 

-fx f n f + ^2{€ k a\a k + V k {fa k + h.c.)} 



where z is the coordination number, ip = (b) the su- 
perfiuid order parameter and V k the hybridization for 
fermions. 

For the phases with a density- wave order, we divide 
the bipartite lattice system into two sublattices [3|| . The 
corresponding Hamiltonian is 



i?GSIAM = 



-zt b (ipabl + p^b a ) + —n b a {n b a - 1) - n b a n b a 



-/-i f n{ + ^{e k al a a ka + V ka (p a a ka + h.c.)} 
k 

+ U b fninl] , 



(G) 



where a = A, B represents the sublattice index (a — A,B 
with a 7^ a). For this Hamiltonian, we perform single- 
site DMFT calculations for each sublattice structure. In 
the following, we fix the parameters, U b = 1.0, zt b = 0.05. 

We calculate the superfiuid order parameter ip and the 
self-energy Y± b f '(u>) self-consistently by employing the nu- 
merical renormalization group (NRG) [36l . l37| as an im- 
purity solver. NRG has the advantage in performing the 
high-accuracy calculation in the low energy region thanks 
to the logarithmic discretization of the conduction band. 
This method has been already extended to include bo- 
son degrees of freedom (24j], which are incorporated in 
the impurity Hamiltonian. This allows us to apply NRG 
with the same accuracy in the low energy region as the 
ordinary fcrmion case. 

We compute several thermodynamic quantities and the 
quasiparticle weight defined by 



Z = 



1 



1 - (ffi£(w)/du)[ 



(7) 



+U bf n b n f , 



(5) 



which is inversely proportional to the effective mass of 
fermions. This quantity represents how strong the cor- 
relation effect is. We also calculate the density of states 
(DOS), p(w), i.e. the single particle excitation spectra 
derived from the imaginary part of the Green's function. 



III. MANY-BODY EFFECTS IN NORMAL 
PHASES OF FERMIONS 

We first discuss many-body effects in the BF mix- 
ture system without a density-wave long-range order of 
fermions. We employ the effective impurity model de- 
scribed by Eq. JS]) for gcncralizcd-DMFT calculations. 



A. Correlation effects in a metallic phase 

In order to figure out possible phases in the BF mix- 
ture system, we calculate several quantities such as the 
fermion filling (n/), the boson filling (n b ), the bosonic su- 
perfiuid order parameter ip and the quasiparticle weight 
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Z as a function of the chemical potential for bosons. In 
particular, the quasiparticle weight Z can be used as a 
measure of correlation effects: for a free particle system 
Z = 1, while for an extremely correlated system Z — > 0. 

The computed results are shown in Fig. [T] Note that 
the origin of the chemical potential Afi b for bosons is de- 
fined so that it gives the fillings (n/) = 1/2, (n b ) = 5/2 
for U bf = 1.0 and [7 b/ = 2.0. At non-integral fillings, the 
boson sector is always in a superfluid phase with finite (p, 
while for (n b ) = 1,2,3,4 (upper panel) and (n b ) = 1,4 
(lower panel), it is in an insulating Mott phase with 
ip = 0. It is seen that there is no renormalization ef- 
fect (Z = 1) without superfluidity (tp = 0). On the other 
hand, in the presence of superfluidity, the quasiparticle 
weight Z of fermions decreases from unity, implying that 
the renormalization of fermions occurs. Figure [T] also 
elucidates that the strength of the renormalization de- 
pends on the magnitude of superfluid order parameter 
(/?; stronger renormalization (smaller Z) occurs for larger 
ip. These results certainly suggest that the many-body 
effects in the fermion sector are induced by the boson de- 
grees of freedom via the boson- fermion interaction U b f . 

We note here that the strong boson-fcrmion interaction 
may possibly induce orderings such as a density wave 
ordering both for the fermion and boson sectors. This 
kind of instability indeed appears around the region near 
Afi B ~ in Fig. [TJ (U bf = 2.0). The corresponding 
data for the physical quantities are lacking there because 
such density-wave ordering is not taken into account in 
the homogeneous DMFT calculations, so that there is 
no convergent solution. We will separately discuss the 
results for the density- wave state in Sec. HVl 

Wc now discuss how the renormalization effects appear 
in the DOS for fermions. In Fig|2l wc show a typical 
profile of the DOS for fermions in a metallic region with 
bosonic superfluidity at U b f = 1.0 and A/i& = 0.0 where 
fillings are (w) = 1/2 and (n ) =5/2. There is a sharp 
peak structure at the Fermi level due to the renormal- 
ization effects with bosonic superfluidity, which is one 
of the characteristic properties of the present BF mix- 
ture system. One may immediately notice that a similar 
peak structure is quite commonly observed in correlated 
electrons in condensed matter physics. Actually, there 
is a close relationship and a crucial difference between 
the present BF system and the electron systems; in both 
cases the strong renormalization of the DOS is caused by 
low-energy collective excitations, but in the former (lat- 
ter) case the collective excitations come from electrons 
themselves (additional boson degrees of freedom). 

Therefore, we naturally expect that some intriguing 
aspects of correlation inherent in the BF systems should 
appear, which are not observed in correlated electron 
systems. Such examples can be indeed found in the 
empty/filled band limit of fermions in our case. We can 
see from FigQ] that the renormalization of fermions oc- 
curs even in the extreme conditions for fermion filling: 
the empty band limit ((tv) ~ 0) or the fillcd-band limit 
((n/) ~ 1). Since these cases provide unique aspects of 



the BF systems beyond ordinary electron systems, in the 
next subsection, we will give more detailed discussions 
on the many-body effects in the two limiting cases. 



U bf =1.0 
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FIG. 1: (Color online) Fermion filling (n^), boson filling (n b ), 
superfluid order parameter ip and the quasiparticle weight Z 
as a function of the chemical potential A/u 6 for bosons for 
fixed U bf = 1.0 (top) and U bf = 2.0 (bottom). 
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FIG. 2: (Color online) DOS for fermions with the boson- 
fermion interaction U b ^ = 1.0 and A^i b = 0.0, where the 
renormalization effects occur in the presence of bosonic su- 
perfluidity. 
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FIG. 3: (Color online) Quasiparticle weight Z as a function 
of the superfluid order parameter ip in the empty band limit. 
The boson-fermion interaction is fixed at U bi = 1.0 (+ points) 
and U bf = 2.0 (x points). 



FIG. 5: (Color online) Quasiparticle weight Z as a function of 
the superfluid order parameter (p at empty band limit. The 
boson-fermion interaction is fixed at U b * = 1.0 (+ points) 
and U bf = 2.0 (x points). 
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FIG. 4: (Color online) DOS as a function of the superfluid 
order parameter ip at U b ^ = 1.0 (left panel) 2.0 (right panel). 
The origin of the DOS is shifted for easy-to-see. 



FIG. 6: (Color online) DOS as a function of the superfluid 
order parameter ip at U b ' = 1.0 (left panel) 2.0 (right panel). 
The origin of the DOS is shifted for easy-to-see. 



B. Correlation effects in empty- and filled- band 
limit 

1. empty-band limit 

We first discuss the relation between the quasiparti- 
cle weight and the bosonic superfluid order parameter 
in the empty-band limit {(n*) ~ 0). The computed re- 
sults are shown in Fig. [3] It is seen that the evolution of 
the superfluid order parameter enhances the renormaliza- 
tion effects. This figure also indicates that the increase 
in the boson-fermion interaction U h f makes the renor- 
malization effects stronger (smaller values of Z). For 
U b f = 2.0, the quasiparticle weight becomes almost zero 
around ip ~ 1.2, which means that the quasiparticle mass 
of fermions becomes extremely heavy. This tendency may 
be related to instability toward the density-wave forma- 
tion, as mentioned above. 

In order to address how the renormalization affects dy- 
namical properties, we calculate the DOS for fermions for 
several choices of ip. The results are plotted in Fig. |4] 
Note that only the particle-addition spectrum has finite 
values (w > 0) because we cannot remove fermions from 
the system in this limiting case. In this figure the DOS 
has a peak near the Fermi level which becomes sharper 
as ip becomes larger. This is consistent with the relation 
between the quasiparticle weight and the bosonic super- 



fluidity mentioned above. 

We note here that this kind of BF mixture in the 
empty-band limit has been already realized in recent 
experiments (T5 , \vj\ , which is sometimes referred to as 
a "polaronic system" [3l|. We hope that the above- 
mentioned renormalization effect could be observed ex- 
perimentally in the near future. 



2. filled-band limit 

We next focus on the opposite extreme case, a filled- 
band limit, where the fermion filling is close to unity 
((nf) ~ 1), while the boson sector is still at non-integer 
filling. In this limit, we can still tune the boson filling 
to control the amplitude of superfluidity. We plot the 
quasiparticle weight Z as a function of the superfluid 
order parameter ip in Fig. [5] Note here again that only 
the particle-removal spectrum has finite values (w < 0) 
because we cannot add fermions into the system in this 
limiting case. 

Figure [5] suggests that the increase in the superfluid 
order parameter enhances the renormalization effects, re- 
sulting in smaller values of Z. Also, we can see that the 
boson-fermion interaction enhances the renormalization 
effects as for the empty band limit. We show typical ex- 
amples of the DOS for several choices of tp with fixed 
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FIG. 7: (Color online) DOS for fermions in the imp 
Hamiltonian with several boson-fermion interactions (L 
0.0, 1.0, 2.0, 3.0) in the presence of the bosonic superfhuuioy. 
A small peak structure on the Fermi surface disappears at 
U — 0.0. Inset shows the enlarged view in the region around 
the Fermi surface. 



interactions in Fig. |6] With increasing ip, the peak be- 
comes sharp in accordance with the corresponding be- 
havior of the quasiparticle weight Z as a function of ip. 
Therefore we conclude that the many-body effects induce 
the peak near the Fermi level in the DOS for fermion sys- 
tems. 

As mentioned above, in the filled-band limit, only the 
particle-removal spectrum can be observed. Experimen- 
tally, this limiting case may be more tractable than the 
empty-band limit if one could use the rf spectroscopy. We 
note that in this case, the renormalization effects should 
be discussed for a "hole" type quasiparticle. 

Before concluding this section, we would like to briefly 
discuss the properties of the impurity Hamiltonian eq. 
(5) in order to exclude the possibility that the peak orig- 
inates from the DMFT iteration process. We use NRG as 
an impurity solver with a constant density of states and a 
constant hybridization in order to extract general proper- 
ties of the impurity Hamiltonian. Filling and interaction 
parameters are (n^) = 1/2, (n) = 5/2, and U b = 1.0. 
Fig. [7] shows the DOS for fermions with several boson- 
fermion interactions in the presence of the bosonic su- 
perfluidity. The DOS has two structures for finite boson- 
fermion interactions (U b f ^ 0). The main part of the 
excitation spectrum comes from the bare hybridization. 
On the other hand, the small peak on the Fermi sur- 
face only appears in the presence of finite boson-fermion 
interactions. We also confirm that the peak structure 
only appears in the presence of the bosonic superfluid- 
ity, which elucidates that the many-body effects occur in 
terms of boson degrees of freedom. These results are con- 
sistent with the DMFT calculations in the lattice system 
and support that the peak structure originates not from 
the DMFT iterations but from the interplay of the boson 
and fermion degrees of freedom. 
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FIG. 8: Schematic phase diagram for the BF Hubbard model 
with fixed fillings {rJ) = 1/2, (n b ) = 5/2. SS means a super- 
solid phase and AMI an alternating Mott insulator. 



IV. MANY-BODY EFFECTS IN A 
SUPERSOLID PHASE 



We now turn to the many-body effects in the phase 
with a density-wave order in the fermion sector. This 
case was already studied theoretically by Titvinidze et al. 
with DMFT [24], [25[ , so we will perform complementary 
calculations to highlight the importance of many-body 
effects which was not addressed in the previous work. 

To treat a density- wave order, we employ Eq.® as 
an effective impurity model, which allows us to treat the 
sublattice symmetry breaking. Here, we present the re- 
sults for fixed fillings (w) = 1/2 and (n b ) = 5/2, where 



we define (n^ b '>) = 1/2 



A,B 



1q (& \ We briefly sum- 
marize the phase diagram with a density-wave order of 
fermions [2J, [25[ . As far as U b f ^ , there always ex- 
ists the difference in the particle density between two 
sublattices, AN b W> = l/2\n a U " > - n h S f) \ ^ 0, signal- 
ing the stability of a density wave state. The density 
wave state is referred to as an alternating Mott insula- 
tor (AMI) for bosons and as a "charge" density wave 
(CDW) for fermions, where "charge" is used following 
the tradition in the solid state physics. Furthermore, 
for AN b m + 1/2 (m is integer), the boson sector 
favors a supcrfluid phase because the commensurability 
condition is not satisfied. Therefore, in the presence of 
bosonic superfluidity, the BF mixture system becomes a 
"supersolid" (SS) phase j25[, which is a main topic in this 
section. 

In order to focus on the essential points, we here fix 
the boson-boson interaction U = 1.0. In this condition, 
we end up with the phase diagram as a function of U b * , 
as shown schematically in Figj8] The physical properties 
in this phase diagram are essentially the same as those 
obtained by Titvinidze et al. [25| for different fillings 
{n') = 1/2 and (n?) = 3/2. Interestingly, they pointed 
out that the DOS for fermions has a peak structure near 
the Fermi level in the SS phase. We demonstrate here 
that this peak is caused by the many-body effects. 

To address the above point, we focus on the supersolid 
SS2 phase where the correlation effects are observed more 
clearly than the SSI phase. We show the DOS in the 
SS2 phase in Fig. |9] It is seen that the DOS has two 
characteristic structures, as expected. In the high energy 
region, there is a hump structure which comes from the 
mean-field type effect of U b f . This exists both in the 
metallic and insulating phases. On the other hand, in 
the low-energy region, there is a sharp peak near the 
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FIG. 9: (Color online) (a) DOS for = 1.8, 1.9, 2.0, 2.2, 2.4. 
(b) DOS enlarged around the Fermi level. 



Fermi level, which only emerges in the supcrsolid phase, 
as already found by Titvinidze et al. for different fillings 
0. 

We now provide evidence that this peak indeed orig- 
inates from the correlation effects due to boson degrees 
of freedom. We show the DOS for different choices of 
the interaction U bf in FiglH With increasing U b f , the 
weight of the peak initially increases and then decreases, 
as seen in the lower panel of FigO This behavior is sim- 
ilar to the one observed between the quasiparticle weight 
and the bosonic superfluidity, so that it is consistent with 
the results obtained in the previous section. In order to 
check the sublatticc-dependence, we show the DOS for 
each sublattice in Fig. [TO] In both cases, the DOS has a 
peak near the Fermi level, suggesting that not only the A- 
but also B-sublattice has the anomalous peak, which was 
not obvious in the previous study [13, HI I- Note that the 
shape of DOS is quite different from each other. In the 
supcrsolid phase, the occupation number of fermions at 
each site is alternating between two sublattices: if AN? 
is close to full ((n^) ~ 1), the B-sublattice remains al- 
most empty ({rig) ~ 0), and vice versa. Note that the 
condition for each sublattice approximately corresponds 
to the empty- and filled band limit discussed in the previ- 
ous section, and therefore the corresponding DOS indeed 
exhibits analogous properties discussed in Figs. 2] and 




FIG. 10: (Color online) DOS for the A and B sublattices at 
U bJ = 2.0. 
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FIG. 11: (Color online) Quasiparticle weight Z as a function 
of superfluid order parameter tp for several values of interac- 
tion 1.6 < U bf < 2.5. Fillings are fixed at (n f ) = 1/2, (n b ) = 
5/2. +, x correspond to the A, B sublattice in Fia llOl 



[6] Therefore, we can say that the emergent peak struc- 
ture in the supersolid phase is a fingerprint of many-body 
effects inherent in the BF mixture system. 

To further confirm our statement, we show the relation 
between the quasiparticle weight Z and the superfluid or- 
der parameter tp in Fig. 1111 The data include the calcu- 
lations performed for several different interactions (1.6 < 
U bf < 2.5) with fixed fillings (n f ) = 1/2, (n b ) = 2/5. It is 
seen that the increase of tp enhances the renormalization 
effects (i.e. smaller Z), whereas there is no correlation 
effect (i.e. Z = 1) in the absence of the bosonic super- 
fluidity tp — 0. This behavior is consistent with that for 
the metallic phase in the previous section. Therefore we 
confirm that the renormalization effects in the supersolid 
phase are induced by the boson degrees of freedom. 



V. CONCLUSION 

We have theoretically investigated a mixture of bosons 
and fermions loaded in an optical lattice using a gen- 
eralized DMFT combined with NRG. We have revealed 



7 



that strong correlation effects emerge in the fermion sec- 
tor in the presence of bosonic superfluidity regardless of 
whether the system is metallic or insulating. This conclu- 
sion has been drawn via a systematic study of the close 
relationship between the renormalization factor and the 
magnitude of superfluidity. 

This kind of boson-driven renormalization effect gives 
rise to the characteristic peak structure in the low-energy 
region of the DOS. The formation of such a many-body 
peak is similar to that expected for the ordinary Fermi 
liquid, but there is a crucial difference between them. 
While in the ordinary Fermi liquid, the renormalization 
is caused by the low-energy excitations of fermions them- 
selves, in the present mixture system, it is driven by 
boson degrees of freedom. Therefore, the latter brings 
about some unique correlation effects inherent in the 
BF mixture systems: for example, the strong renormal- 
ization appears even in the empty /filled band limit of 
fermions. Since the empty band limit was already re- 
alized in recent experiments, we hope that many body 
effects proposed here will be observed in the near future. 

We have confirmed the appearance of the peak struc- 
ture in the impurity model and concluded that the many- 
body effects occur in the presence of the bosonic super- 
fluidity even in the impurity Hamiltonian. This implies 
that the peak structure originates from the boson degrees 
of freedom, not from DMFT iterations. In the preced- 
ing work [25] an instability towards the phase separation 
was pointed out as a possible mechanism for the peak- 
formation. Although in the present calculation, we have 



not encountered such a tendency, more detailed analyses 
should be necessary to figure out the relationship between 
our scenario and the previous one. 

Although the calculation in this paper has been done at 
T = 0, we expect that such an anomalous peak structure 
in the DOS can be observed with the rf-spectroscopy ex- 
periments j28H3l| at sufficiently low temperatures where 
fermions and bosons are both in the quantum degenerate 
regime. The emergence of such a peak structure should 
be a fingerprint of the many-body effects inherent in the 
BF mixture. 
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